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Abstract. Non-abelian energy loss in quark gluon plasma is shown to lead to 
novel hadron ratio suppression patterns in ultrarelativistic nuclear collisions. We 
apply recent (GLV) estimates for the gluon radiative energy loss, which increases 
linearly with the jet energy up to E < 20 GeV and depends quadratically on the 
nuclear radius, R. The K~*~/7r + ratio is found to be most sensitive to the initial 
density of the plasma. 


Energy loss of high energy quark and gluon jets penetrating dense matter 
produced in ultrarelativistic heavy ion collisions leads to jet quenching and thus probes 
the quark-gluon plasma formed in those reactions [[jj [|]. The non-abelian radiative 
energy loss, A E(E,L), suppresses the moderate 2 — 3 GeV < pr < 10 — 15 GeV 
distributions of hadrons in a way that can also influence the jet fragmentation pattern 
into different flavor hadrons. This is because quark and gluon jets suffer different 
energy losses proportional to their color Casimir factors (4/3, 3). First estimates |j], 
suggested that A E ss 1 — 2 GeV(L/fm) would depend linearly on the plasma thickness, 
L , as in abelian electrodynamics. In BDMS || non-abelian (radiated gluon final state 
interaction) effects were shown to lead to a quadratic dependence on L with a much 
larger magnitude of A E. The analysis of jet-quenching in Pb + Pb at yfs = 20 AGeV, 
however, indicated a negligible energy loss at that energy ||]. In GLV JJ j7j finite 
kinematic constraints were found to reduce greatly the energy loss at moderate jet 
energies and we obtained A E(E, L) ~ E ■ (L/ 6 fm) 2 , Here we apply the GLV energy 
loss to estimate more quantitatively the effect of jet quenching on the hadronic ratios 
at moderate p± < 10 GeV at the Relativistic Heavy Ion Collider (RHIC) energy, 
y/s = 130 AGeV. We illustrate the jet quenching effects for a generic plasma with 
an average screening scale p = 0.5 GeV, strong coupling a s = 0.3, and an average 
gluon mean free path X g = 1 fm. We introduce the opacity as the mean number of jet 
scatterings, n = L/X , which is assumed to be a finite number. 

An advantage of looking into particle ratios is that uncertainties in the absolute 
normalization due to acceptance tend to cancel. The disadvantage is of course that 
high p± particle identification is increasingly difficult. As we have shown in Ref. ||, 
the kinematically suppressed GLV energy loss turns out to depend approximately 
linearly on the jet energy, E. This linear dependence, A E oc E, leads to only a very 
weakly p± dependent suppression of the transverse momentum distributions for the 
kinematic range accessible experimentally at RHIC. We focus on the p± dependence 
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of the particle ratios in order to help pin down the jet quenching mechanism. Our 
primary candidate is the measurable K + /n + ratio. 

In order to investigate the influence of the GLV energy-dependent radiative energy 
loss on hadron production, we apply a perturbative QCD (pQCD) based description 
of Au + Au collisions, including energy loss prior to hadronization. First, we check 
that the applied pQCD description reproduces data on pion and kaon production in 
p+p collision. Our results are based on a leading order (LO) pQCD analysis. Detailed 
discussion of the formalism is published elsewhere &E3- 

Our pQCD calculations incorporate the parton transverse momentum (“intrinsic 
hr") via a Gaussian transverse momentum distribution g(kr) (characterized by the 
width (/cf,)) H [lO], [yj with the usual convolution of the parton distribution functions 
(PDF) f a / p , partonic cross sections and fragmentation functions (FF) D h / C , as: 
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Here we use LO PDF from the MRST98 parameterization [|l2| and a LO FF |13| |. 
The applied scales are Q = pd^ and Q = pr/'JZc. 

Utilizing available high transverse-momentum (2 < pr < 10 GeV) p + p data on 
pion and kaon production at 19 < yfs < 63 GeV we can determine the best fitting 
energy-dependent (kj) parameter (see Ref. §0 for further details on (fcf,)). One can 
conclude about the K + /tt + ratio that, while agreement is reasonable for pr > 4 GeV, 
there is a systematic discrepancy at smaller transverse momenta. High statistics p + p 
data at RHIC will be essential to establish an accurate baseline to which A + A must 
be compared. Now we turn to the calculation for Au + Au collision at RHIC energy, 
\/s = 130 AGeV. As a first approximation, let us consider slab geometry, neglecting 
radial dependence. We include the isospin asymmetry and the nuclear modification 
(shadowing) into the nuclear PDF [ p~4| ]. 

The value of the (kj) of the transverse component of the PDF will be increased 
by multiscattering effects in A + A collisions. Two limiting cases were investigated 
with (i) a large number of rescatterings [|ll| and (ii) a small number of rescatterings 
(“saturated Cronin effect”) |j| [ffj. Here we consider the influence of jet-quenching 
on the hadron spectra, leaving the inclusion of the multiscattering effect for future 
work. In this simplified calculation, jet quenching reduces the energy of the jet before 
fragmentation. We concentrate on y cm = 0, where the jet transverse momentum 
before fragmentation is shifted by the energy loss, p*(L/ A) = p c — A E(E,L). This 
shifts the z c parameter in the integrand to z* = z c /( 1 — A E/p c ). The applied scale 
in the FF is similarly modified, Q = pt/2z*, while for the elementary hard reaction 
the scale remains Q = p c /2. 

With these approximations the invariant cross section of hadron production in 
central A + A collision is given (in a somewhat abbreviated notation) by 
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where B = 2n / 0 &max b db Taa( b) with b max = 4.7 fm for the 10 % most central Au+Au 
collision and Taa( b) is the thickness function. The factor z* r Jz c appears because of 
the in-medium modification of the fragmentation function |lq| . Thus, the invariant 
cross section (|j) will depend on the average opacity or collision number, h = L/X g . 

We note that the GLV energy-dependent energy loss leads to a rather structureless 
downward shift of the single inclusive yields of all hadrons because AE/E is 
approximately constant. This is in contrast to the results in JnJ, where an energy 
dependent fractional energy loss, A E/E = (0.5 GeV / E)(L/fm), was used. 

More generally, the quenching factor depends also on fluctuations of energy loss 
that occurs because a single collision dominates the energy loss. We are presently 
testing the effect of fluctuations using the normalized f(SE) oc d(SE — SE 0 )/SE 2 
distribution of energy loss from |(| and SE 0 fixed by ( SE) = AEqlv{E,L). 
Preliminary results indicate, that fluctuations simply renormalize the effective 
opacity that is needed to obtain a fixed quenching factor. The reason that this 
simplification occurs in spite of the large effect of fluctuations is the approximate 
energy independence of the GLV fractional energy loss in the relevant jet energy range. 
Since other effects, such as expansion, also renormalize the value of the effective static 
opacity, we continue here to illustrate the effects of jet quenching simply varying L/\ g 
as discussed above. 

To understand the origin of the different dependence of the tt and I\ quenching 
factors on L/\ g , we consider the gluon contribution to the hadron yield for charged 
pions and kaons. For simplicity, we illustrate the gluon contribution to the non- 
quenched ratios only. Fig. 1. displays (top panel) the relative contribution of gluon 
fragmentation to the production of pions and kaons as a function of pr- It can 
be seen that gluons dominate pion production while (valence) quarks dominate I\ + 
production, especially at lower pt values. Due to the dependence of jet energy loss 
on the color Casimir factor, pions will experience stronger suppression than kaons for 
the same value of h = L/X. There is a pr dependence of this effect because the quark 
and gluon fragmentation functions differ. While 7r + and tx~ show an almost identical 
relation between quark and gluon fragmentation contribution, there is a difference 
between K + and A'” in accordance with the different production mechanism of K + 
and K~ [jn|. Our main result is the prediction that there could be a factor of 2-3 
enhancement in the K + /tt + ratio in the pt ~ 3 — 6 GeV range for effective opacity 
h = 4. The recently reported suppression of the 7r° spectra in Ref. 0 corresponds 
to an effective static opacity in the range h ~ 3 — 4 [Q . 

In summary, we applied the GLV parton energy loss [pj, Q] in LO pQCD to pion 
and kaon production in nuclear collisions. The energy dependence of A Eqlv leads 
to a weak pr dependent but flavor dependent quenching factor that increases with 
the effective opacity of the plasma. The enhancement of the K + /ir + ratio in the 
2 < p± < 10 GeV range was found to be the most sensitive to the gluon opacity at 
RHIC energies. 

Acknowledgments 

This work was supported by the U.S. DOE under DE-FG02-86ER40251, DO-AC03- 
76SF00098 and DE-FG02-93ER40764, by the U.S. NSF under INT-0000211, by the 
BCPL in Bergen and OTKA No. T029158, T032796. 


K + /tt + Probes of Jet Quenching in AA Collisions 


4 


o 

Z5 

O' 


-h 


cr> 



O 


cn 


CH 




p T (GeV) 

FIG. 1. Top panel: Relative contribution of gluon fragmentation to charged pion and kaon 
production at y/s = 130 AGeV without jet quenching. Bottom panel: /iv + ratio in 

Au + Au —> h + X collision at yfs = 130 AGeV without jet-quenching (n = 0) and with it, 
n = 2,3,4. 
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